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Recent theoretical insights into the possibility of non-Abelian phases in o — 2/3 fractional quantum Hall 
states revived the interest in the numerical phase diagram of the problem. We investigate the effect of various 
kinds of two-body interlayer couplings on the (330) bilayer state and exactly solve the Hamiltonian for up to 
14 electrons on sphere and torus geometries. We consider interlayer tunneling, short-ranged repulsive/attractive 
pseudopotential interactions and Coulomb repulsion. We find a 6-fold ground-state degeneracy on the torus 
when the interlayer hollow-core interaction is dominant. To identify the topological nature of this phase we 
measure the orbital-cut entanglement spectrum, quasihole counting, topological entanglement entropy, and 
wave-function overlap. Comparing the numerical results to the theoretical predictions, we interpret this 6-fold 
ground-state degeneracy phase to be the non-Abelian bilayer Fibonacci state. 
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Introduction. There is a surging interest in the search for 
non-Abelian (NA) anyons in topological states of matter [1- 
8]. NA anyons register protected quantum memory and the 
stored information can be processed through braiding the 
world line of NA anyons. This comprises the basic notion 
of building fault-tolerant quantum computers. Two famous 
examples are Ising and Fibonacci anyons which have been 
conjectured to emerge in i/ = 5/2 and i/ = 12/5 fractional 
quantum Hall (FQH) states, respectively [9-11]. In partic¬ 
ular Fibonacci anyons are of interest for universal quantum 
computation. However, there is no experimental observation 
of these exotic quasiparticles in FQH systems to date. An¬ 
other approach for realizing NA anyons is through coupling 
Abelian states via various types of interactions to drive phase 
transition to NA topological phases [12-17]. In a realistic ex¬ 
perimental situation, some of the interactions considered in 
theoretical investigations are indeed relevant [18-25]. Pursu¬ 
ing this direction further can potentially lead to a novel venue 
for realizing NA anyons. 

In this Rapid Communication, we focus on the two- 
component FQH system at ^ = 2/3, whose parent state con¬ 
sists of two decoupled z/ = 1/3 Laughlin states. The two com¬ 
ponents may describe physical spins [18], spatial (physical 
layer) degrees of freedom, e.g., in the wide quantum well real¬ 
ization of the 2/3 state [19-22], or valley indices in graphene¬ 
like systems [23-25]. For convenience, we refer to all these 
various realizations of two-component systems as the “2/3 bi- 
layef’ state. Early studies of the 2/3 bilayer state focused on 
various possible Abelian phases [26] and the experimental ob¬ 
servation of phase transitions [18, 19]. In particular McDon¬ 
ald and Haldane numerically established the phase diagram 
of the system in the presence of interlayer tunneling and the 
bare Coulomb interaction projected to the lowest Landau level 
(LLL) [27] and confirmed the possibility of various Abelian 
phases. Recent analytical insights into the possibility of NA 
phases in the 2/3 bilayer state [3, 28-33] brought this seem¬ 
ingly closed problem into new focus. However, the results 
of Ref. [27] imply that microscopic realizations of these NA 


phases require perturbation to the model Hamiltonian consist¬ 
ing of interlayer tunneling and the bare Coulomb interaction 
projected to the LLL. In fact such perturbations may already 
exist in experimental systems. 

Therefore, the recent proposals [3, 28-33] beg for revisiting 
the 2/3 bilayer problem numerically considering the modified 
interlayer Coulomb interaction and interlayer tunneling. Any 
two-body interaction projected to the LLL can be expanded in 
terms of the so-called Haldane pseudopotentials, Vm [34, 35], 
which are projector operators of two electrons to the state with 
relative angular momentum m. Modifying the Coulomb in¬ 
teraction amounts to changing the coefficient of the pseudo¬ 
potentials. Here, we modify the interlayer Coulomb interac¬ 
tion by changing its dominant components, namely and 
ymter known as hollow-core interaction). We solve the 
resulting Hamiltonian for up to 14 electrons using the exact- 
diagonalization (ED) method and utilize a variety of numer¬ 
ical measurements to establish the nature of different topo¬ 
logical phases achieved through varying coupling parameters. 
We obtain a wide range of parameters where enhancing 
and suppressing Vg"'" components of the Coulomb interaction 
drives phase transition to a NA state with Fibonacci anyons. 
Model. We start with two v — 1/i Laughlin states that are 
coupled through a number of distinct interlayer interactions. 
Let us now consider a torus with Ly dimensions. In the 
Landau gauge A = B{0, x), single-particle states in the LLL 
are labeled by their momentum along the y direction quan¬ 
tized as ky = 27rnjLy. The parent state, the (330) Halperin 
state [36], is the exact eigenstate of the following Vi Hal¬ 
dane’s pseudopotential 
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where cr =t,4' denotes the layer index, c„ annihilates 
an electron with ky = 27rn/Lj, momentum, Vr]} = 
(r^ — s^) 2 " , and n = 2Tt/Ly. 

We investigate the effects of interlayer coupling by consid- 
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ering both tunneling and interlayer two-body interaction of 
Haldane’s Vq and Vi pseudo-potentials. The tunneling Hamil¬ 
tonian is 


= -t± E 


5 


( 2 ) 


where = (c„ c„ ^). The two-body interactions are 

Vr^ = Uo Vr^°Jci^,^A+s,-.Cr,,-aCn+r+s,., 0) 
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In order to address the stability of the phases that we ob¬ 
tain with this model we also consider the effect of interlayer 
Coulomb interaction in the limit of negligible interlayer sepa¬ 
ration perturbed by the above terms. 

Candidate phases. Below, we consider two categories of 
possible states ?A u = 2/3: Abelian and non-Abelian, each 
with three members. Table I summarizes some of their defin¬ 
ing properties, (i) The (330) state, which is simply two decou¬ 
pled :/ = 1/3 Laughlin states. This phase has 9-fold degener¬ 
ate ground states on the torus, (ii) 1/3: The particle-hole(PH) 
conjugate of the 1 / = 1/3 Laughlin state [37, 38]. (iii) Layer 
singlet state: The nonchiral (112) Halperin state [36-39]. (iv) 
The bilayer Fibonacci phase [3], which is a non-Abelian state 
described by SU ( 3)2 Chern-Simons field theory and contains 
6 distinct quasiparticles, three of which are non-Abelian with 
d = quantum dimension, (v) The interlayer (1. L.) 

Pfaffian state [32], which has 9 distinct quasiparticles and 
its wave-function can be written as the product of the ( 221 ) 
Halperin state and a Pfaffian wave-function [32]. (vi) The 
fermionic (M = 1) Read-Rezayi (RR) state [28, 29, 40]. 
This non-Abelian state with 15 distinct quasiparticles is de¬ 
scribed by Z 4 X U{1)q edge conformal field theory. A unique 
feature of this state is the emergence of a non-Abelian quasi¬ 
hole excitation with — e /6 fractional charge. 

Numerical results. Now we numerically search for possible 
phases by exact diagonalization of our model Hamiltonian. To 
mitigate the limitations of finite-size system study, we make 
measurements of various characteristics of topological phases 
using both toroidal and spherical geometries. The number of 
electrons A/ and total number of quantum fluxes are re¬ 
lated through the relation on the torus and on the 

sphere — S, where S is the shift of a FQH state 

on the sphere [41]. Good quantum numbers that can be used 
to label energy states are the total center-of-mass momentum 
K on the torus and angular momentum on the sphere. In 
the absence of the interlayer tunneling, the total pseudospin 
Sz = i(-^T ~ ^ good quantum number, where 

is the electron number in the upper (lower) layer. 

Below we present the resulting phase diagram and the mea¬ 
surements characterizing the phases. 


Phase c GSD 5 V 

( 330 ) 2 9 3 3 

173 0 3 0 Vs 

(112) 0 3 1 Vs 

Fibonacci 14/5 6 3 (5 -f %/5) 

I.L. Pfaffian 5/2 9 3 ^12 

^4 RR 2 15 3 6 


OES 

1,2,5,10,20,--- 

1,1,2,3,5,--- 

1,1,3,6,13,--- 
1,2,6,13,--- 
1,1, 3, 6,13,-•• 


Table I. Total central charge c, ground-state degeneracy (GSD), 
the sphere shift 5 = — N^, total quantum dimension 

D, and orbital-cut entanglement spectrum (OES) counting for I = 
0,1, 2, 3,4, • • ■ (I is the angular momenta relative to the ground 
state) for various topological phases at v — 2/3. For NA states, 
the OES counting generally depends on the number of electrons 
and the pseudospin in region A. Here we only show the case 
when j 2 (N^/4) is an integer and 5'7 = 0 for the Fibonacci 
and I. L. Pfaffian states (Z^ RR state). Although the countings of the 
bilayer Fibonacci state with = 0 (mod 2) and Z 4 RR state with 
= 0 (mod 4) are equal for Z < 4, they are different for other 
values of N^. The counting of the (112) state is the convolution of 
the counting of two counter-propagating U{1) modes. 



Figure 1. (Color online) The phase diagram based on GSD of the 
Hamiltonian H = -f Z7r“V(,"‘“ -f -f Uf The num¬ 

bers in the parentheses indicate GSD. (a)-(c) The GSD for Aj, = 8 
particles for (a) = 0, (b) L]"'" = 0, and (c) = 0. (d) 

The GSD as a function of 17]"’" for = 8 and 10 electrons with 
jjmter — =0. Eq,Ei, E 2 are the lowest three energy levels in 

the K < A^/3 sectors. One can see that with the increase of 17]"'“, 
E 2 — El ^ El — Eq, implying a total 6-fold degeneracy. 


(a) Ground-state degeneracy. The GSD on the toms is 
significant as it equals the number of distinct anyon excita¬ 
tions [42]. Figure 1 shows the phase diagram of the interaction 
H = -F -F + Ht exhibiting several 

phases with 3-fold (purple), 6 -fold (green), and 9-fold (red) 
GSD phases. Of these the 9-fold GSD region in Figs. 1(a) and 
1 (b) is smoothly connected to the parent state, so it is most 
likely the (330) phase. For the 3-fold GSD region, there are 
two Abelian candidates (see Table I) which can be differen¬ 
tiated by their different shifts on sphere, S. The most tan- 
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Figure 2. (Color online) Here we consider the interaction F/^ouLmb + 
"Vf" wi^^ = -0.4and(7i"“ = 0.6 

on the square torus [43]. (a) The energy spectrum of = 6,8,10, 
and 12 electrons in the = 0 sector, (b) The finite-size scaling of 
the energy gap and the ground-state splitting in the = 0 sector, 
(c) The PES of = 8 electrons; the counting below the gap (green 
line) is 1332. (d) The PES of = 10 electrons; the counting below 
the gap (green line) is 4875. 


talizing is the 6 -fold GSD region, for which the NA bilayer 
Fibonacci state is the only candidate to the best of our knowl¬ 
edge. Figures l(a)-l(c) show that dominant is required 
for the 6 -fold GSD phase. Nevertheless, this phase is stable 
against subdominant tunneling and terms. In the rest 
of this Rapid Communication, we will focus on the = 0 
limit where we can take advantage of as a good quantum 
number. Figure 1(d) shows that this degeneracy is stable upon 
increasing the system size. 

To further address the stability of the 6 -fold GSD phase 
against higher order pseudopotentials which must exist in any 

realistic setting, we study the H = iT'Sufomb + -f^couTomb + 
[jmterymter [/mteryinter Hamiltonian with < 0. We Con¬ 
sider this particular form, because the component of the 
bare Coulomb interaction is larger than its component. 
Intriguingly, Fig. 2(a) shows that this Hamiltonian can sta¬ 
bilize the 6 -fold GSD phase with a robust degeneracy and a 
large gap. Moreover, the finite-size scaling of the energy gap 
and the ground-state splitting [Fig. 2(b)] shows that this phase 
is probably gapped in the thermodynamic limit. The finite- 
size scaling of the energy gap on the sphere also supports this 
conclusion [Fig. 3(c)]. 

(b) Edge state counting and entanglement spectrum. The 
orbital-cut entanglement spectrum (OES) [44] encodes the de¬ 
generacy of the edge excitations with angular momentum I 
relative to the ground state [39]. We compute the OES on the 
sphere which leads to a single boundary. Eor a bilayer EQH 
system without interlayer tunneling, each OES level on the 
sphere can be labeled by the total electron number N^, the 
total angular momentum L^, and the total pseudospin 5^ in 
the subregion A. 


Table I lists the expected counting for the candidate phases 
that can be obtained using the thin torus patterns and gener¬ 
alized exclusion rules [45, 46]. Care has to be given to the 
fact that OES counting generally depends on N^. Eor exam¬ 
ple, in the case of Abelian phases, the OES is independent of 
N^, while for the bilayer Eibonacci state it depends on the 
parity of N^. In the case of the Z 4 RR state, the counting 
depends on modulo 4. In general, it can be shown that 
the OES counting of a fc—electron clustered topological phase 
described by SU{n)^. Chern-Simons theory depends on 
modulo k. Eor the bilayer Fibonacci state, the counting for 
= 2M + 1 is 1, 2, 5, 10, • • • . For the RR state, the 
counting is 1, 2, 5, 10, • • • (1, 2, 6 ,11, • • •) for = AM — 1 
(N^ = AM — 2) for large values of N^. For finite systems, 
the counting is slightly different. For instance, when = 6 
the expected counting for the bilayer Fibonacci (Z 4 RR) state 
is 1,1, 3, 6 ,12, • • • (1,2, 5, 7, • • •). For = 7, the expected 
counting of the RR state becomes 1, 2,4, 7, • • •. 

For the region with 9-fold GSD we obtained OES counting 
consistent with the Abelian (330) state: 1,2,5,10, 20, • • • in¬ 
dependently of N^. On the other hand, the measured OES 
counting for the 6 -fold GSD phase is different for the odd and 
even values of and distinct from any Abelian state. The 
OES counting of the modified interlayer Coulomb interaction 
with 6 -fold GSD is presented in Figs. 3(a) and 3(b). The OES 
counting we obtain for the bilayer Fibonacci state dictates the 
following generalized Pauli exclusion rules: (a) In every three 
consecutive orbitals, there are at most two electrons, (b) Each 
configuration of two electrons with distance more than two or¬ 
bitals such as ( 1001 ) or ( 10001 ) is doubly degenerate after as¬ 
signing spin indices and must be counted twice. It was shown 
that these rules give rise to 6 -fold GSD on the torus, and re¬ 
sult in three NA anyon excitations with d = quantum 

dimension, hence the Eibonacci anyons [3]. 

The fact that the counting for ( = 1 is unambiguously 1 
for = 6 rules out the I. L. Pfaffian and Z^ RR states (for 
both of these states at ( = 1 , we expect the counting to be 2 ). 
Moreover, the fact that the measured OES counting depends 
on the parity of indicates that the 6 -fold GSD phase is a 
paired state with k = 2. This is consistent with the SU{3)2 
Chern-Simons description of the bilayer Eibonacci state. 

(c) quasihole counting. The degeneracy of the ground-state 
in the presence of quasihole excitations, usually referred to 
as quasihole counting, is a telling property of a topological 
phase [48, 49]. One way to create quasiholes is by adding 
additional fluxes to the system. Eor instance, in the bilayer 
Fibonacci phase, each additional flux would create two quasi¬ 
holes. Another way is to use the particle-cut entanglement 
spectrum (PES) which is achieved by tracing out some elec¬ 
trons in the density matrix associated with the ground-state 
wave function [50]. Each PES level is labeled by the total 
electron number N^, momentum Kj^, and pseudospin Sf in 
the subsystem A. In Eigs. 2(c) and 2(d), we observe a clear 
gap in the PES. The counting below gap is different from the 
expected counting for the Abelian as well as I. L. Paffian and 
Z 4 RR states [46], while it is consistent with the generalized 
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Pauli exclusion rules of the bilayer Fibonacci state described 
in the previous section. 

(d) Topological entanglement entropy. The entangle¬ 
ment entropy of a two-dimensional gapped state follows the 
5'^ (L) = aL — 7 relation, where L is the perimeter of sub- 
region A. The subleading constant 7 is related to the total 
quantum dimension of the topological phase T), through the 

7 = log I? [51, 52] relation. The precise numerical eval¬ 
uation of the topological entanglement entropy (TEE) 7 for 
small systems is usually challenging [47, 53-55]. However, 
qualitative aspects of our results in Eig. 3(d) suggest that the 
6 -fold GSD phase is a NA state. The numerically estimated 7 
for the 6 -fold GSD state is larger than the estimated 7 for the 
unperturbed (330) state and far from the values of 7 for the 
remaining Abelian phases whose 7 ’s are expected to be half 
the value of the (330) state (see Table I). Although this argu¬ 
ment suggests the resulting state is NA, it does not rule out the 
possibility of other NA states. 

(e) Wave-function overlap. We have also computed the 
squared overlap of the ground state of the 6 -fold GSD re¬ 
gion on the sphere with several candidate states in Table I. The 
overlap with the (112), 1/3, and Z 4 RR state is always neg¬ 
ligible. The overlap with the (330) state is only 0.12 for 12 
electrons. The overlap with the I. L. Pfaffian state is 0.80 for 

8 electrons, which is not quite low, but it decays very quickly 
to 0.73 and 0.66 when the system increases to 10 and 12 elec¬ 
trons. 

To gain analytical insight about the numerical results, we 
also studied the thin torus limit of the coupled (330) system 
[46]. We show that at least four of the topological phases 
conjectured in Table I are plausible. Eurther, we obtain the 
thin torus patterns of the ground states [56-59]. We then 
use the resulting patterns to extract several key pieces of in¬ 
formation about the nature of the underlying topological or¬ 
der, e.g., GSD, fusion rules, and quantum dimensions of the 
anyons [3, 58]. In particular we confirmed that the numeri¬ 
cally obtained ground states of the 6 -fold GSD phase in the 
thin torus limit are consistent with the thin torus patterns of 
the bilayer Eibonacci state [46]. 

Discussion. In summary, we reinvestigated the phase dia¬ 
gram of the coupled (330) bilayer state and considered the ef¬ 
fects of the modified interlayer Coulomb interaction by vary¬ 
ing its and components. We have verified the 

result of McDonald and Haldane [27] indicating that inter¬ 
layer tunneling and bare Coulomb interaction projected to the 
ELL can only yield Abelian phases. However, we found that 
when the pseudopotential components decay smoothly and 
both Pi = and 52 = Vfincrease, the re¬ 

sulting ground-state exhibits many properties of the bilayer Ei¬ 
bonacci phase. Thus we have provided numerical evidence for 
the topological quantum phase transition between the (330) 
state and a NA phase. Moreover the robustness of the bilayer 
Eibonacci phase which contains Eibonacci anyons which are 
capable of performing universal quantum computation is tan¬ 
talizing. 

Our result raises the question of how the 6 -fold GSD re- 


Figure 3. (Color online) For (a)-(c) we consider the modi¬ 
fied Coulomb interaction with HlfSomb + -f 

jymterymter jymter ^ _q ^ fainter ^ q g 

(b) show the OES for = 14 electrons in the Nf — 6, — 0 

and Nf = 7, Sf = —1/2 sectors. The levels consistent with the ex¬ 
clusion rule proposed in the text are colored by red. (c) The finite-size 
scaling of the energy gap in the = 0 sector, (d) The extrapolation 
of the measured entanglement entropy Sa{L) to the thermodynamic 
limit, S'tl, for different boundary lengths L, using the method in 
Ref. [47] for the unperturbed (330) state and the 6-fold GSD phase 
in the large-Vf“' limit. 


gion of our model can be realized experimentally. Our calcu¬ 
lations indicate that when and p 2 ratios are large enough, 
the Fibonacci state becomes energetically favorable. If the 
interlayer distance is small, so pi « 1. One 

strategy to get larger p 2 is to go to higher Landau levels with 
2/3 filling [60, 61]. It is known that in higher Landau levels, 
the Vg component of the bare Coulomb interaction decreases 
with respect to ELL while other components do not change 
too much. Although in the second Landau level (2LL) for two 
dimensional electron gas (2DEG) systems, p 2 « 0.68 is sig¬ 
nificantly larger than its value 0.5 in the ELL, unfortunately 
we do not have smooth pseudopotentials (the V 2 component 
is larger than the V 4 component). However, in the 2LL for 
graphene-like systems, p 2 « 0.70 grows and pseudopoten¬ 
tials decay smoothly. In fact, it can be shown through ED 
that even the bare Coulomb interlayer interaction projected to 
the 2 LL in graphene yields a 6 -fold degenerate ground-state, 
and considering slight finite interlayer distance stabilizes that 
phase even further [62]. Therefore, we expect that graphene 
is a promising platform for the Fibonacci state. 

Note added. Recently we became aware of related works in 
Refs. [63-65]. 
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SUPPLEMENTAL MATERIALS 

In this supplemental material, we study the thin torus limit of our model Hamiltonian for 2/3 bilayer fractional quantum Hall 
(FQH) system analytically in section A. We show that when hollow-core interaction dominates, we expect phase transition in 
the bilayer Fibonacci state. In section B, we present a method based on Ref. [44] to derive the orbital-cut entanglement spectrum 
(OES) counting for the bilayer Fibonacci state with an even number of electrons in subregion A. A similar calculation would 
give the counting for other situations and remaining phases. In section C, we compare the particle-cut entanglement spectrum 
(PES) counting of the ground state in the 6 -fold degeneracy region with various candidate states at ly = 2/3. Einally in section 
D, we show the squared wave-function overlap of the ground state in the 6 -fold degeneracy region with various candidate states. 


A. Thin torus limit of Hamiltonian 


In this section we study the (330) EQH state in the thin torus limit (Ly/L^ <C 1) [56-58] and focus on four regimes, in each 
of which one interlayer coupling dominates over others and we can ignore other interactions. In the Landau gauge two copies 
of the Laughlin state atv= 1/3 filling [called as the (330) state] is the exact eigenstate of the following Hamiltonian (usually 
referred to as Vi Haldane’s pseudo-potential [34]) in the second quantization representation 


V]"' 


= EE 

n,r,s 




“'n+r+s,<7 5 


(5) 


where a =t 5 i denotes the layer index, c„ g. annihilates an electron with ky = ‘I'KnjLy momentum, Vr}} = 

3 2 T’^ I - 2 

(r^ — s^) e~'^ 2 , and k = 2tt jLy. In this supplementary material, we consider three types of interlayer couplings: 

(a) Interlayer tunneling with 


coupling Hamiltonian, (b) Interlayer Haldane’s Vq pseudo-potential. 


V inter_ 

0 — 


t/o ^ V}] 

(T]n,r,s 

u„ •£ 


(0) t t 

■5 ^n+r,(T^n+s,——(7^n+r+s,cr5 




n,r,s 


(7) 


where VrJ = 2 ^ and s„ ^ ~ Cn.^Cm.t) '^^e destruction operator for the (layer-)singlet state 

between sites n and m (Note that n and m indices are related to the momentum of electrons along y direction. In the Landau 
gauge, they also dictate the average position of electrons along x axis). In this letter, we consider both repulsive {Uq > 0) and 
attractive ((7q < 0) cases, (c) Hollow-core interaction, i.e., interlayer Haldane’s Vi pseudo-potential: 


yinter^ 


Ui 


( 1 )J J 


‘-'n+r,<7^n+s, ——cr^n+r+s,cr 


(7;n,r,s 


t/i E 


J 


n+s,n+r^n,n+'r+S5 


( 8 ) 


where f„ ^ ^ destruction operator for a (layer-)triplet state between sites n and m. In the thin 

torus limit, the parent Hamiltonian as well as interlayer interaction greatly simplify and evaluate to : 


Vf 




^ ^ ^ ^ ^r,0 '^n+r,a't 
n r—1,2 


(9) 


^oE 




0.0 “f ^,0 ^n,n-\-r^n,n-\-r 


( 10 ) 


yintei- 




^1 / . '' 1,0 ^n,n+A'n,n+t^ 


n 


( 11 ) 





7 


Phase 

GSD 

Thin torus patterns 

(330) 

9 

[ti, 0 , 0 ], [t, 4-1 0 ], [4-, ti 0 ]-(-translations 

1/3 

3 

[—>■, —>, 0 ] - 1 - translations 

( 112 ) 

3 

['^’Ts 0 ]+translations 

Bilayer Fibonacci 

6 

[t 4 -i 0 , 0 ], [^”‘ 4 )^’'*^, 0 ]-(-translations 

Interlayer Pfaffian 

9 

[t-[, 0 , 0 ], [^’'‘y2^’'^,0], ['^’“’'‘-^^-’“’'^[-Htranslations 

Z 4 Read-Rezayi 

15 

[->■, ->■, ->■, -)■, 0 ,0], [->■, ->■, ->•, 0 , -)■, 0 ],[-(-, ->■, 0 , -)■, -)■, 0 ]-(-translations 

3-Pfaffian state of charge 2e bosons 

9 

[ 74-1 0,0, t4-i 0 , 0 ], [t 4 -i 0 , til 0 , 0 , 0 ]-(-translations 


Table II. Ground-state degeneracy and the thin torus patterns associated with some of the candidate phases at t/ = 2/3. The number of different 
patterns for each phase equals the corresponding GSD. In the thin torus limit, ground-state is represented by repetition of either patterns. The 
resulting patterns encode the spatial density of electrons. Electrons with t (/) label live on top (bottom) layer. —>■= (t + i) represents 
resonating electrons. 


where ^ = Cn,crC„ „ is the electron number operator at site n. V”'” is already diagonal in the occupation basis and can 
be easily solved. for spin up sector contains three ground-states: |t, 0 ,0, f, 0, 0, f, 0,0, • • •), |0, f, 0,0, f, 0,0, f, 0, • • •), 
|0,0, t, 0, 0, t, 0,0, tj • ■ ■ Similarly the spin down component of V™'” can be solved and the ground-state of full intra-layer 
Hamiltonian is simply the 9 different possible tensor products of the ground-states of the parts, (see Table II). 

Below we study the four extreme limits: 

(i) |fx| ^ \Uq \ ,Ui. In this limit, we must first diagonalize , and then consider the effect of intra-layer repulsion. In order 
to minimize "Hj, we need to project spin along +x direction (if f is a real positive number), i.e. electrons must occupy the 
symmetric state |—)')j = ^ (|t)^ -b |j,)j). Then there are two regimes: 

(a) is much larger than then we obtain the following ground-state patterns: 

[—)■, —>■, 0 ] = |—)•, —>■, 0 , —>■, —>■, 0 , 0 , •••), 

[ 0 , —>■, — >■] = jo, —>■, —>■, 0 , 0 , , 

[—>■, 0 , —>■] = 0 , —>■, —>■, 0 , 0 , —>■,•••). ( 12 ) 


By consulting Table II, we conclude that for very extreme interlayer tunneling we expect the particle-hole conjugate of the 1/3 
Laughlin state 1/3 state. 

(b) t± is comparable or less than then using second order perturbation we obtain the following ground-state patterns [3]: 


(t, 4 - + 4 -, t) , 0 ] + two other translations, 

V2 


(13) 


where as before [a,b,c] denotes the building block (unit cell) of the ground-state. The ground-states can be obtained by repeating 
these unit cells all over ID lattice associated with occupation number basis. Table II suggests it is (112) Halperin state [36]. 

(ii) Ui 3 > |f_L|, \Uo\- In this case, we again must deal with Since, it consists of triplet counting operators, the ground-state 
must be spin singlet between near neighbors. Doubly occupied sites are also spin singlet. Therefore, the thin torus patterns of 
the ground-state patterns are: 


[ 200 ] = [n, 0 , 0 ], 


[ 110 ] ^ [ 


1 

71 


(t, 4 - — 4 -) t) j 0 ] + translations. 


(14) 


In this way, we obtain six degenerate ground-states. The thin torus patterns tells us that every three consecutive orbitals contain 
exactly three orbital in the ground-state. According to table II it represents the bilayer Fibonacci state which contains Fibonacci 
NA anyons. 

(iii) Uq ^ Ui'. Considering the onsite term in because of strong repulsion between two layers, ground-state cannot 
contain doubly occupied sites. This leaves us with six low energy patterns. However, taking the nearest neighbor interaction, i.e. 
UqVi qs\ favors triplet formation. Therefore, we expect the following three ground-states: 


[-7 (f, j') ^ 0] -I- translations. 

v2 


(15) 


According to table II the above thin torus patterns are signatures of (112) Abelian state. 

(iv) —Uq ^ j/j.!: Ui'. When interlayer coupling is attractive and strong, we expect electrons on the same form a bound-state. 
These bound-states can be imagined as charge — 2e bosons. Assuming these electrons pairs are localized and tightly bound, we 
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can project the remaining terms in the Hamiltonian to the this low energy subspace and obtain an effective interaction, between 
bosons. Doing so, we obtain the following effective Hamiltonian: 

-> E E (^r[o -\U°\ Kp ) nUA (16) 

n r—1,2 

where = blbi is the density of bosons (electron pairs). Now consider the following two groups of patterns: 

(а) : [ti, 0,0, ti, 0,0] + translations, 

(б) : [ti, 0, ti, 0,0,0] + translations. (17) 

These two groups are not exactly degenerate with respect to the effective Hamiltonian Eq. (16) in the thin torus limit. However, 
by considering appropriate three-body interactions to the effective Hamiltonian Eq. (16), the two groups above can become 
degenerate and evolve together. In that case, the resulting Hamiltonian describes nine degenerate ground-state patterns. These 
patterns are indeed root configurations of the ground-state wave function. Comparing the above nine patterns with those in Table 
II, we conclude that the above phase with additional three-body interactions is described by the 3-Pfaffian state. 


B. Edge state counting using the thin torus patterns 


In this section we closely follow the method developed in Ref [44] to compute the edge state or entanglement counting of 
the bilayer Eibonacci state. The counting of other states can be found similarly. Eor this end, let us consider the half infinite 
cylinder geometry with a boundary at x = 0. In the thin cylinder limit, one of the ground-states of the bilayer Eibonacci state is 
the following pattern (this is the unique ground-state pattern on sphere geometry): 


|GS)i 


n, 0 , 0 , n, 0 , 0 , n, o, o, n, o, o, o, o, o, o, o • 


(18) 


On the left side of a; = 0, the (fj,, 0,0) pattern is repeated all the way to —oo, while there is no electron on its right side. The 
total momentum (or total angular momentum on sphere geometry) of an arbitrary state is defined as: 

i = E = E (19) 

iez a 

Let us define I as the relative momentum of an arbitrary state with respect to the ground-state in Eq. (18), i.e. I = L — Cqs- 
example, I = 0 for the ground-state itself. Now we want to study excited states with relative angular momentum 1. There are 
two categories of excitations: 

(a) Edge excitations: These are low energy excitations which respect the generalized Pauli exclusion rules for the bilayer 
Eibonacci state. 

(b) Bulk excitations: Any state which does not respect the generalized Pauli exclusion rules. 

The Pauli exclusion rule for the bilayer Eibonacci state is: (i) In every three consecutive orbitals, there should be at most two 
electrons, (ii) Electrons with distance one or two (e.g. 110 or 101 configurations) from singlet pairs and must be counted once. 
The formation of 110 singlets has priority unless the formation of 101 singlets can increase the total number of singlets, (iii) 
Electrons with distance more than two (e.g. 1001 or 10001) can form both singlet and triplet states and thus must be counted 
twice. 

In the following, we use the generalized Pauli exclusion rule mentioned above to find the counting of the edge excitations for 
( = 1,2, 3,4. 

Example 1. ^ = 1: There is only one admissible edge excitation (the states in this section need to be normalized): 

|( = 1, edge)i = • • • n, 0,0, ti, 0,0, U, 0,0, (f, i - i,t), 0,0,0,0 • • • (20) 

Example 2. I = 2: There are three admissible edge excitations: 

\l = 2, edge), = • • • n, 0,0, fi, 0,0, fi, 0,0,0, U, 0,0,0,0, • • • 

\l = 2, edge)2 = • • • n, 0,0, fi, 0,0,0, (f, i t), 0, (t, i t), 0,0,0,0, • • • 

\l = 2, edge)3 = • • • n, 0,0, fi, 0,0, fi, 0,0, (f, 0, i - i, 0, f), 0,0,0 • • • , 

Example 3. ( = 3: There are five admissible edge excitations: 

|/ = 3, edge), = • • • n, 0,0, U, 0,0, U, 0,0,0, (f, ; - i, t) , 0,0,0, • • • 

|/ = 3, edge )2 = • • • n, 0,0, U, 0,0, (f, i t) , 0,0, U, 0,0,0,0, • • • 


( 21 ) 
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(a) (b) 



Figure 4. (Color online) The OES of N^. — 14 electrons in the = 6 , S'^ = 0 and = 7, S'^ = —1/2 sectors in the limit of 
^^intcr ^ _ g countings that suffer from finite-size effect are labeled by green numbers. 


|/ = 3, edge)3 = • • • n, 0,0, U, 0,0, (f, i - i, t), 0, (f, 04 - i, 0, f), 0,0,0, • • • 

|/ = 3, edge)4 = • • • 4 ,0,0, (f,; - i, t) , 0, (t, i - 14) , 0, (t, i - 14) , 0,0,0,0, 

|( = 3, edge )5 = • • • 4,0,0,4,0,0,4,0,0, (f, 0,04 - i, 0,04), 0,0, • • • ■ 

|/ = 3, bulk)e = • • • 4, 0,0, 4, 0,0, 4, 0,0, (t, 0,04 + i, 0,04), 0,0, • • • • 

Example 4. ( = 4: There are ten admissible edge excitations; 

|( = 4, edge), = • • • 4, 0,0, 4, 0,0, 4, 0,0,0,0, 4, 0,0,0, • • • 

|( = 4, edge)2 = • • • 4,0,0,4,0,0,0,4,0,0,4,0,0,0,0, • • • 

|( = 4, edge )3 = • • • 4, 0,0, 4, 0,0, 4, 0,0,0, (f, 04 - i, 04) , 0,0, • • • 

|( = 4, edge )4 = • • • 4, 0, 0 , (t, i - 14) , 0, (t, i - 14) , 0, 0 ,4, 0 , 0,0,0, • • • 

|( = 4, edge )5 = • • • 4,0,0, (f, i - 14), 0, (t, i - 14), 0, (t, 04 - 0, f), 0,0,0, • • • 

|/ = 4, edge )6 = • • • 4, 0,0, 4, 0,0, (f,; - 4) , 0,0, (t4 - 4) , • • • 

|/ = 4 ,edge )7 = • • 4i, 0,04i, 0,044,0,044, • • • 

|/ = 4, edge )8 = • • 4i, 0,0,4,0,0, (f, 0,; - i, 0, t), (t, 0,; - i, 0, t), 0,0,0, • • • 

|/ = 4, edge)g = • • 4i, 0,0,4,0,0, (f,; - i, t), 0, (t, 0,0,; - i, 0,0, t), 0,0, • • • 

|/ = 4, bulk) 10 = • • 4i, 0,0,4,0,0, (t, i - i, t), 0, (t, 0,0, i + i, 0,0, t), 0,0, • • • 

|/ = 4, edge)ii = • • 4i, 0,0, 4, 0,0, 4, 0,0, (f, 0,0,0, i - i, 0,0,0, f) , 0, • • • 

|( = 4, bulk) 12 = • • 4i, 0, 0 ,4, 0 , 0 ,4,0, 0 , (t, 0 , 0 , 0 , i o, o, o, t), o, • • • 

|( = 4, edge)i 3 = • • • (t4 - 14) , 0, (t4 - i, t) , 0, (t, i - 14) , 0, (t, i - i, t) , 0,0,0,0, • • • (24) 

Note that |( = 4, edge ),3 can happen only for > 8 , and thus the counting for = 6 and I = 4 is 12 rather than 13. 
Similarly, we can find the counting for higher angular momenta in the = 0 (mod 2) and = 0 sector. The OES counting 
for other values of Nf modulo 2 and Sf can also be obtained with possible slight changes of the exclusion rule. In the limit of 
^mter ^ ymtra jjig counting is Very clean (Fig. 4). 


( 22 ) 

(23) 


Phase 

4 = 8 ,4^ = 4 

4 = 10,iVf = 4 

4 = 12 , 4^=6 

(330) 

784 

3025 

48400 

( 112 ) 

784 

3025 

48400 

Interlayer Pfaffian 

1125 

6085 

104491 

Z 4 Read-Rezayi 

210 

715 

7506 

6 -fold degenerate state 

505 

2242 

32840 


Table III. The PES counting of the ground state in the 6 -fold degeneracy region as well as various candidate states on the sphere for different 
system sizes in the = 0 sector. 


C. Quasihole counting using the PES 

In this section we show the PES counting of the ground state in the 6 -fold degeneracy region as well as various candidate 
states on the sphere (Table III). We can see that the counting of the ground state in the 6 -fold degeneracy region is very different 
from all candidate states that we have checked, thus ruling their possibilities out. 
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Phase 

IVe =6 

00 

II 

N, = 10 

IVe = 12 

(330) 

0.44 

0.28 

0.18 

0.12 

Interlayer Pfaffian 

0.89 

0.80 

0.73 

0.66 


Table IV. The squared wave function overlap of the ground state in the 6 -fold degeneracy region with various candidate states on the sphere. 
We use the Hamiltonian iTcoufomb + -^rauiomb + -I- with = —0.4 and t/J"'" = 0.6 to get the ground state in the 6 -fold 

degeneracy region. We use H — Vl"“ + -I- T-Lt with = 2 to get the 1 /3 state. We antisymmetrize the (330) state to get the Z 4 state. 

Other states are generated from their parent Hamiltonians. The overlap with the (112), 1/3 and Z 4 RR state is always negligible. 


D. Wave function overlap 

We have computed the squared wave function overlap of the ground state in the 6 -fold degeneracy region with various candi¬ 
date states on the sphere (Table IV). The overlap with the (112), 1/3 and Z 4 RR state is always negligible. From the aspect of 
the overlap, we can see that the main competing phase is the inter-layer Pfaffian. However, the ground-state degeneracy on the 
torus as well as OES and PES can rule it out. 










